International Journal of Mechanical and Production 
Engineering Research and Development (UMPERD) 

ISSN(P): 2249-6890; ISSN(E): 2249-8001 
Vol. 5, Issue 6, Dec 2015, 81-116 
© TJPRC Pvt. Ltd. 

MATLAB PROGRAM FOR TRADITIONAL DESIGN AND 
OPTIMIZATION OF 30 KW SCIM MOTOR 

RAVI PRASAD 1 2 & K PRAHLADA RAO 3 

1 Retired Scientific Officer (F), Bhabha Atomic Research Centre, Mumbai, Maharashtra, India 
2 Former In-charge Head, Department of Nuclear Energy, Pandit Deendayal 
Petroleum University, Gandhinagar, Gujarat, India 
3 Professor, Department of Mechanical Engineering, JNTUA College of 
Engineering, Anantapur, Andhra Pradesh, India 

ABSTRACT 

Various steps in development of a MATLAB program for the traditional design of 30 KW SCIM motor have 
been detailed. The whole process of traditional design of motor design is divided into five parts and again each part 
into various steps. The MATLAB program listed performs all the activities associated with five parts and their steps of 
the process. Relevant comments are added into the MATLAB program at suitable locations to identify the calculations 
associated. Parameters of design of motor and the equivalent variables used in the MATLAB program of steps of 
design process under each part are listed in the relevant tables. The list of variables used and equivalent MATLAB 
variables completely bring clarity to the steps followed in the traditional design of the motor. 
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I INTRODUCTION 

Several methods of reduction of losses have been identified and reported in report [3] by the same authors 
of this report. As part of the thermal design improvement of Squirrel cage induction motor (SCIM) motors which 
forms the main focus of the thesis work of author 1, this MATLAB program is developed to study the effects of 
reduction of losses due to the techniques adopted and detailed in the above report. 

In this report step by step traditional procedural design based evaluation of motor along with optimisation 
both for round conductors and bar conductors of a 30 KW SCIM motor has been presented. The MATLAB 
program listing is provided as the annexure. In the text book [1] two MATLAB programs, one for the evaluation 
of motor parameters and the other for optimisation of parameters for a 30 KW motor are explained. This 
MATLAB program is an improvement compared to the two MATLAB programs combined of text book [1] from 
the following aspects. In the present report both of these have been combined into one single program, with a 
value of variable OPTIM being given a value equal to 1 (motor optimisation ) or Option 0 (No optimisation). 

II MATLAB PROGRAM FOR TRADITIONAL DESIGN AND OPTIMISATION OF MOTOR 

Various steps in the overall design of motor using MATLAB program have been explained with tables 
and equations and the corresponding MATAB expressions have been indicated with suitable comment statements. 
For a 30 KW motor many manufacturers still prefer to have round conductors; at the same time usage of bar 
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conductors also exists. So the provision of designing with both types of conductors is made available in the program. When 
round conductors are used; it is a common practice to use tapered slots and straight teeth. The slot area in the program 
corresponds to the mean area of tapered slots assuming a fixed taper for all the slots. Mean height to mean width ratio of 
the slot is a parameter of variation of optimisation. And slot fill factor has been the other parameter of optimisation. Report 
[2] gives the sequential steps provided in the traditional design of a motor but the traditional design is done using hand 
calculations which do not cover the scope for automation of design of motor. The values of detailed calculations have been 
provided for bar conductors (cond_type=2) and the output values are provided for rod conductors (cond_type=l) with 
optim variable value equal to 0. Optimisation results of the program for rod conductors and bar conductors are given in 
tables XXXXI and XXXXII. Optimisation of slots with different stator and rotor slot profiles and air gap lengths for 
number of laminations (n_lams) also has been made available in this program for a given value of n_lams. In the present 
program provision for n_lams has been made for a value of 4 but the values have been given same and n_lams is given 
equal to 1 . 

Table 1 gives the listing of details of variables corresponding to motor parameters used in the mathematical 
expressions and equivalents of MATLAB program. 


Table 1 


Part I 

MATLAB 

Used for 

Variable used 

Variable 

in Equations 

SKW 

KW power of motor 


SBav 

Average flux density 

Bav 

SKWa 

Kilo watt power of motor 


ac/m (sq) 

Electrical power density 

q 

Vph 

Volts per phase 

v Dh 

f 

frequency 

f 

FI 

Flux tp 

<p 

ki 

Iron factor 

... 

CO 

Output coefficient 

Co 

SPF4P 

Power factor values at the 
given rated power 


SEF4P 

Efficiency values at the 
given rated power 


KW 

Kilo watt of motor 

Q 

V 

Voltage 

V 

Iph 

Current per phase 

I D h 

PP 

Pole Pitch 

ttD 
T ~ P 

LbyPP 

Length/Poe Pitch 

L/t 

P 

Number of poles 

P 

V 

Peripheral velocity 

V 


Table 2 gives the values and units of the parameters corresponding to name plate details of a motor. 


Table 2: Name Plate Details of 3 Phase 4 Pole SCIM Motor 


Parameter 

MATLAB 

Variable 

Variable 

Units 


Rated power 

KW 

Q 

[KW] 

30 

Phase voltage 

V 

v Dhl 

V 

690 

Input phase 
current 

Iph 

Iphl 

A 

27.07 


Impact Factor (JCC): 5.6934 


NAAS Rating: 2.45 
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Table 2: Contd., 

Synchronous 

speed 

Ns 

N s 

rpm 

1500 

Supply frequency 

f 

f 

[Hz] 

50 

Targeted power 
factor: 

pf 

COS((> 

-- 

0.84 

Targeted 

efficiency: 

eff 

n 

-- 

0.927 

Lamination 

thickness 

DELT 

A 

mm 

0.5 

Shaft diameter 

riry 

Try 

mm 

110 

Sots/pole/phase 

spp 

q s /q, 

— 

4 

Winding factor 

Kw 

Kpdl 

— 

0.955 


A. Limitations in Design 


Table 3: Selection Criterion for Physical Dimensions of Motor 



Range of ( L/r ) 

Will Yield/ Result 

a) 

1.0 to 1.1 

Good overall design 

b) 

1.0 to 1.3 

Good power factor 

c) 

1.4 to 1.6 

Good efficiency 

d) 

1.5 to 2.0 

Minimum overall cost 


The materials used for the machine imposes a limitation in design. The limitations stem from saturation of iron, 
current density in conductors, temperature, insulation, mechanical properties, efficiency, power factor etc. 

Ill STEP 2 OF PART I - DESIGN CRITERION FOR BASIC DIMENSIONS OF STATOR 

Main physical dimensions of motor are obtained to achieve the design criteria which are obtained based on ratio 
of gross length and pole pitch Ratio (L/r) as explained below: 


Table 4: Physical Dimensions of Motor 


Step 2of Part I Physical Dimensions of Motor 


MATLAB 

variable 

Variables used in 
Equations 


Diameter of Stator 

D 

D 

0.213 

Gross Length 

L 

L 

0.207 

D 2 L 

DsqL 

D 2 L 

0.0112 

Net Iron Length 

Li 

L, 

0.190 


Table 5: Average Performance of Class 2 Motor 


Speed 

Performance at rated output 



Torque 

Nm 

Current 

n 

Pf 

1465 

195 

58 

91.6 

0.82 

n 

Lor DL 
stator Pu 
I S /In 

Break down 
torque T s /T n 

T max /T n 


91.6 

6.6 

2.4 

2.8 



Efficiency Class = 2 T s = Starting torque I s = Starting current 


III STEP 3 OF PART I BEARING AND FRAME 

The frame sizes are generally standard and are made of either A1 alloy or Cl. A1 alloys are better conductors of 
heat where as cast iron material has higher strength. Motors made of better tolerances and lower values of air gap have 
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better heat dissipation from rotor to stator and then to environment. 

Table 6: Bearing and Frame Details of Motors 


Frame 

Details 

Bearings 

Radial Ball Bearing 
Details 


DE 

Side 

NDE 

Side 

Bore 

OD 

3Width 

200 L 

6312 

6312 

60 

130 

31 


IV. STEP 4 OF PART I - MAGNETIC AND ELECTRIC LOADS SELECTION 

Iron losses largely depend upon air gap flux density 

Table 7: Advantages of Higher Value of B AV 


Size of the 

Cost of the 

Overload 

Machine 

Machine 

Capacity 

Reduces 

Decreases 

Increases 


Table 8: Selection Criterion for Physical Dimensions of Motor 


ac/m 


The Value of Electric Loading Varies 
between 5000 and 45000 i.e. 5000 < ac < 
45000 


Interpolate specific electric loading ac/m (ampere 
conductors per meter of air gap circumference) from 
Table 8. The value of ac/m determines how heavily the 
electric material is utilized. High ac/m means less 
electric material but higher electric losses. 


Table 9: Dependence of Iron Losses on Air Gap Flux Density 


B av 

B Lies between 0.35 and 0.53 i.e 0.35 < B < 
0.53 and the Suitable Values Can be 
Interpolated from the KW Value of the 
Motor Selected. 

Bmax 

The Maximum Flux Density Commonly 
Employed is between 2.2T to 2.4T. 

Select B g from experience (limited by losses in the teeth 
and magnetizing current). B g depicts how heavily the 
magnetic core material is utilized. High B g means less 
magnetic material but higher magnetic losses. Select 
magnetic material also based on frequency. 


V. STEP 5 OF PART I - VALUES OF INTERPOLATION 


Tables X, XI and XII give the values of B av and q, efficiency and power factor values selected from Design data 
hand book for the given motor power. 


Table 10: Values of B A v and ac/M for 30 KW Motor 


KW 

1 

2 

5 

10 

B av (T) 

0.35 

0.38 

0.42 

0.48 

ac/m (q) 

16000 

19000 

23000 

25000 

KW 

20 

50 

100 

500 

Bav (T) 

0.48 

0.50 

0.51 

0.53 

ac/m (q) 

26000 

29000 

31000 

33000 


Impact Factor (JCC): 5.6934 


NAAS Rating: 2.45 
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Table 11: Power Factors vs. Motor Power 


SKWa 

5 

10 

20 

50 

100 

200 

500 

SPF4P 

0.85 

0.86 

0.88 

0.90 

0.91 

0.92 

0.93 


Table 12: Values of Efficiency for Different Motor Powers 


SKWa 

5 

10 

20 

50 

100 

200 

500 

SEFF4P 

0.85 

0.87 

0.88 

0.90 

0.91 

0.93 

0.94 


Table 13: Values of Current Densities -Vs Rotor Diameters 


Diameter 
in Meters 

0.1 

0.15 

0.2 

0.3 

0.4 

0.5 

0.75 

1 

J A/mm 2 

4 

3.8 

3.6 

3.5 

3.5 

3.5 

3.5 

3.5 


Table 14: B-H Curve for 0.5mm LOHYS Material 


B 

.1 

.2 

.3 

.4 

.5 

.6 


H 

50 

65 

70 

80 

90 

100 


B 

.7 

.8 

.9 

1 

1.1 

1.2 

1.3 

H 

110 

120 

150 

180 

220 

295 

400 

B 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 

2 

H 

580 

1000 

2400 

5000 

8900 

15000 

24000 


Carters’ coefficients are available for the ratio of the slot width and the air gap vs carter coefficient and are tabled 
as below. 


Table 15: Carters Coefficients 


Ratio 

0 

1 

2 

3 

4 

5 

6 

For open slots 

0 

.18 

.33 

.45 

.53 

.6 

.66 

For closed slots 

0 

.14 

.27 

.37 

.44 

.5 

.54 

Ratio 

7 

8 

9 

10 

11 

12 


For open slots 

.71 

.75 

.79 

.82 

.86 

.89 

For closed slots 

.58 

.62 

.65 

.68 

.69 

.7 


Table 16: Flux Density vs WPKG (Power Loss of Watts Per Kg) 


B 

0.8 

1.2 

1.6 

2.0 

2.4 

WpKg 

7 

15 

24 

34 

50 


VI. OUTPUT COEFFICIENT (C 0 ) ESTIMATION — GOVERNING EQUATIONS 

Output of a SCIM motor is also equal to 

Q = 3 V phl x I phl x Cos (q>) x q x 10' 3 KW (1) 

Or equation (2) can also be written as 

Q = 3(4.44 x K pdl x f x (p m x N ph i x I ph i x Cos (tp) x q x 10 " 3 (2) 

(Since V phl = 4.44 x K pdl x f x <p m x N phl ) and 
(p m = B av xtxL = B av x j-xL 

Total No of conductors on Stator = 3x2 N ph i = 6 N ph i 

Total Ampere Conductors on Stator = 6 N ph | I ph | which is known as total electric loading 
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Specific electric loading: It is defined as electric loading per meter of periphery, denoted by ac 


ac 


ZNphi'phi __ „ , nD ac 

or IS-nt, 1 In 


phi ‘phi ~ 6 

By substituting the values of f, <p m & N Phl I Phl in equation (2) 

Q = 3 x 4.44 x 0.955 x( x ( B ave x^XL)x x Cos tp x q x 10' 3 KW 

Q = (17.4 x 10" 5 B ave ac Cos tp x q ) D 2 Ln s 
Or Q = C 0 D 2 Ln s 

Where C 0 = Output Co-efficient = ( 17.4 x 10" 5 B ave ac Cos tp x q) 


(3) 

(4) 


Output equation relationship between electrical rating and physical dimensions (Quantities) or the D~L equation 
relating output power (Q), synchronous speed (ns) and rotor volume through an output coefficient (C 0 ) is as given in Eq. 
(1) below: 


C n =- 


kw s n s D z L. 

The VA (Q) rating is proportional to 

• The average gap flux density B; 

• The armature surface current density ac; 

• The volume of the rotating member, D 2 L product and the speed of rotation, n s . 


(5) 


Table 17: Electric and Magnetic Parameters of Motors 



MATLAB 

variable 

Variable used 
in Equations 

Units 


Rating 

SKW 

Q 

(KW) 

30 

Volts 

Vph 

v ph 

Volts 

398 

Step 5 Interpolated 

Values from Tables I to III 

Average flux 
density 

SBav 

Bav 

tesla 

0.484 

ac/m (q) 

q 

q 

amp 

cond/m 

26802 

Efficiency 

eff 

q 

- 

0.892 

Power factor 

Pf 

Cos(tp) 

- 

0.885 

Step 6 Motor Parameters 

Output 

coefficient 

Co 

Co 

Ws/m 3 

107.61 

Synchronous 

speed 

ns 

n s 

rps 

25.00 

Peripheral 

Speed 

V 

V 

(m/s) 

16.73 

MAXIMUM PERMISSIBLE IS 30 M/SEC 


Because the contains air gap flux density (B g ) and surface current density (Ki) in Eq. (2), where a is the magnetic 
saturation factor, K w i is the winding factor. 


Impact Factor (JCC): 5.6934 


NAAS Rating: 2.45 
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D 2 L 


equation has no relationships connecting air gap quantities with the flux density and current densities existing 


in the motor’ s interior. 


Table 18: Effect of Slots on Design 


Less slots 

• 

Less cost 

• 

Less space lost due to insulation and slot opening 

More slots 

• 

Smaller leakage inductance and larger breakdown torque 

• 

Small MMF harmonics 

• 

Better cooling 


VII. DESIGN OF STATOR (PART II) 

Step 1 of Part II 

In this step all the stator dimensions corresponding to the stator poles are estimated. 


Table 19: Variables Corresponding to Design of 

Stator and Equivalent MATLAB variables 


Step 1 of Part II 

sp 

Slot Pitch 

yss 

Tph 

Turns per phase 

N ph 

CDSW 

Current Density 

j 

j 

Current per phase 

Iph 

inSW 

Width wise insulation 

Win 

hstrip 

Height of strip 

hstrip 

HW 

Height of wedge 

h w 

HL 

Height of Lip 

Hi 

Zph 

Conductors per phase 

Z D h 

Zsh 

No of conductors height wise in 
slot 

Zsh 

insH 

Height wise insulation 

h ms 

Zs 

Conductors per slot 

Z s 

Step 2 of Part II 

Rt 

L max 

Max flux density of the tooth 

Btmax 

Rph 

Resistance/phase 

Rph 

L mt 

Length of mean turn 

L m t 

p 

A cus 

Copper loss 

p 

A cu,s 

Step 3 of Part II 

He 

Height of core 

h c 

FIc 

Flux in core 

9c 

Do 

Core outer diameter 

D 0 

Ac 

Area of core 

A c 

Step 4 of Part II 

PitpKg 

Iron loss in tooth/kg 

P,,’ 

Pit 

Iron loss in tooth 

Pi, 

Wt 

Weight of tooth 

w t 

PicpKg 

Iron loss in core/kg 

Pic’ 

Pic 

Iron loss in core 

P,C 

Wc 

Weight of core 

W c 
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VIII. DESIGN OF SLOT SIZE 

Fill factors of 0.4 and 0.5 are used for stator slots when round conductors are packed in slots, and outer diameter 
of slots is calculated where as this could be close to 1.0 for rectangular cross section conductors. Number of conductors per 
slot is found such that flux density is within permissible limits of the material used. Semi-closed and trapezoidal (with 
rounded corners) stator and rotor slots with rectangular teeth are designed for this motor. 

The influence of number of stator slots is as given below table. 

IX. STEP 2 OF PART II- STATOR SLOT GEOMETRY 

There are two types of slots namely: a) Partially open (semi) and b) Closed Slot 

For small motors, taper on the tooth or slot is significant and tapered slots (and parallel sided teeth) are used. This 
gives maximum area of slot for given tooth flux density. Round wires of small gauge are used since they are easy to wind 
and do not mind the taper of the slot. 


Table 20: Details of Conductor Arrangement in Stator Slot 


Slots 

Si 

SI 

48 

Slot Pitch (mm) 

yss 

7tD/Si 

13.941 

Slot pitch provided is close but less than ideal range of 15 to 
20 mm. The slot pitches provided are acceptable for 
partially (semi) closed slots (Tapered slots and parallel 
teeth) 

Conductors/Slot 

N cl 

Zs 

21 

SLOT-PITCH(MM) (PERMISSIBLE: 18T025MM) 

Turns/Ph 

Tph 

N Dh 

120 

No. of conductors /phase 

Zph 

Zph 

240 

Flux/Pole(Wb) 

Fipp 

(p/p 

0.01566 

Phase Current (A) 

Iph 

Iph 

31.78 

STEP 2 OF PART 

II 

Bare Strip (w*t)mm 

Hst,Tst 

w st ,t st 

5.50X1.55 

Width to Thickness Ratio 

hstbytst 

W st ,/t s , 

3.55 

(Permissible:2.5to3.5) 

Area of CS cond(mm 2 ) 

As 

A s 

8.244 

Current density! A/mm 2 ) 

CDSW 

j 

4.000 

No of Strips (width*depth 
-wise) 

Zsw 

Zsw 

3X5 

Slot-Width (mm) 

Wsl 

W sl 

9.5 

Slot-Height(mm)Total 

hsl 

hsi 

43.0 


X. CROSS SECTIONAL AREA OF STATOR CORE 

The cross section area of stator core A cs is given by 

A cs = (Om/2)/ B cs Here, Bcs is flux density in stator core, which is taken about 1.2 Tesla. 

Stator teeth width Value selected should be between 6.4 mm and 25.42 

Ratio of slot width to slot pitch Value selected should be between 0.4 and 0.6 - Ref of Lipo 

Depth of stator core d cs = A cs /L (6) 


Impact Factor (JCC): 5.6934 


NAAS Rating: 2.45 




MATLAB Program for Traditional Design and Optimization of 30 KW SCIM Motor 


89 


Outer dia of stator. Do = (D+2.dss + 2.dcs) (7) 


Table 21: Core Details of Stator 


O.D. of Stator 

Do 

D 0 

mm 

330.0 

Copper c/s in slots 

Acs 

A cs 

mm 2 

8.366 

No. of conductors 
/slot 

Zs 

z s 

- 

18 


XI. STATOR SLOT DIMENSIONS 

Area of one stator slot A ss = Z ss /space factor (Space factor lies between 0.25 and 0.4) 

For parallel sided slot of stator, minimum tooth width of stator, is given by- 

W ts mu, = (P. ®m)/(l .7 L. Ss) (Here, maximum allowable flux density in stator tooth is 1.7 T) 

XII. SLOT SPACE FILL FACTOR FOR OBTAINING SLOT DIMENSIONS 


High voltage machines have lower space factors due to large thickness of insulation. The slot should not be too 
wide to give a thin tooth. The width of the slot should be so adjusted such that the mean flux density in the tooth lies 
between 1.3 to 1.7 Tesla. Width of tooth should not be too large as it results in narrow and deep slots. Deeper slots lead to 
larger value of leakage reactance. In general the ratio of slot depth to slot width should be between 3 and 6. 


Per phase stator current, I ph i = 


Q X io 3 
3V phl Cosip Tj 


F 


cl 


tph.1 


( 8 ) 

(9) 


F cl = area of stator conductor 


Current density § lvalue is obtained from the interpolated value of plot from data available for Stator diameter vs 
current density values 


XIII. STEP 3 OF PART II - FLUX DENSITY ESTIMATION IN STATOR TEETH 


For estimation of flux density in stator teeth various details of stator are estimated at a distance of 1/3 of tooth 
height from narrow end of stator tooth. 


Maximum flux density in stator tooth should not exceed 1.8T; otherwise iron losses and magnetizing current will 
be abnormally high. (So if flux density >1.8T, change slot dimensions) 

^rd 

Mean flux density in the stator tooth is calculated at —of tooth height from the narrow end of the stator tooth. 

1 rd 

Diameter of stator at —of tooth height from narrow end 

D| ht =D + i/i s x2 (10) 


Slot pitch at —of tooth height from narrow end 




T l. = ■ 
sg~ h t 


(11) 
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Width of the tooth at -^-of tooth height from narrow end 


hi, = t i -b q 

t~h t sg-h t s 


( 12 ) 


Area of one stator tooth at —of tooth height from narrow end 
= b t i h x K t l (Where I = k l=Actual iron length) 

Area of all the stator teeth under one pole 

= Area of one tooth x (No of teeth per pole or ^ ) = b t x Kj 



n (o+| h s x 2) 


- b. 


x Ki l X - 

1 p 


(13) 


(14) 


So mean flux density in teeth 


B 


<p 1 


t~h t 


(15) 


Table 22: Slot Details at 1/3 Height from Narrow End 


Slot Details at 1/3 Height from End 

Diameter Equation (10) 

D13 

D 1 

5 h t 

241.6667 

Slot Pitch Equation (11) 

spl3 

T 1 u 

15.8170 

Width Equation (12) 

Wtl3 

b 1 
H h t 

6.2670 

Mean Flux Density 
Equation (15) 

B13 

B 1 

1.0933 

Tooth Flux Density 

Max(T) 

Btma 

X 

Bt maX 

1.6400 


XIV. STEP 4 OF PART II - COPPER LOSSES IN STATOR WINDING 

Mean length of turn of stator winding (L mtl ) 

L mtl = 2 L + 2.3 t +0.24 (16) 

Resistance of stator winding per phase (R Ph i) 

R phl = 0.021 xl0' 6 ^x N phl (17) 

(Only bar conductors are used) Total copper loss in stator winding= 3 I ph 2 R ph i (18) 

Table 23: Copper Losses in Stator Winding 


Length of Mean-Turn (m) 

Lmt 

L mt l 

1.039 

Resistance/Ph (ohm) 

Rph 

Rphi 

0.3175 

depth of Stator Core(mm) 

He 

he 

30.50 

Outer Diameter of Stator Core(mm) 

Do 

D 0 

360.0 

Stator Copper Loss(W) 

Wcus 

l 2 Rphi 

961.9 


P l t = Iron loss in teeth = p’ i>t * 7860 * volume of iron teeth 
P i c = Iron loss in core/yoke = p\ c * 7860 * volume of iron yoke 


Iron loss in yoke = = P i y * Weight of iron in yoke 


Impact Factor (JCC): 5.6934 


NAAS Rating: 2.45 
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Table 24: Iron Losses in Stator 


Interpolated Values of W/KG 

Stator Teeth 

pitpkg 

P\ 

r l,t 

24.91 

Core 

picpkg 

P’ 

r i,c 

18.32 

Weight of Stator (KG) 

Teeth 

wt 

W t 

19.21 

Core 

wc 

w c 

46.83 

Total 

wst 

W st 

66.04 

Iron 

Losses in Stator 

Teeth (W) 

Pit 

Pi,. 

478.6 

Core (W) 

Pic 

p„ 

858.1 

Total 

Pisto 

P i.s.to 

1336.7 


XV. STEP 5 OF PART II -CHECK FOR FLUX DENSITY IN STATOR TEETH AND YOKE 

Table 31 is used for checking the limits. Stator laminations are made of Silicon steel in all the three motors and 
the stator core loss coefficient corresponding to lamination steel is interpolated as per table 15 corresponding to 50 Hz 
frequency. Flux through stator yoke is equal to half of the flux through flux per pole 

Thus B y = — — — and flux density in yoke should be in the range of 1.3 to 1.5 T 

2 l h y Ki 

Corresponding to flux densities in core and tooth of stator B y and B i find out values (p’ ijt and p' ly ) iron loss per 

Z 3 n ' t 

Kg from the graph or interpolate the values corresponding to table 

XVI. STEP 1 OF PART III - ROTOR DESIGN 

Air gap provided is fixed in the MATLAB program but the program is structured such that there is a provision for 
varying the air gap value and optimizing the slot parameters for achieving further optimization of the motor. 

XVII. STEP 2 OF PART III NUMBER OF ROTOR SLOTS 


Table 25: Dimensional Details of Rotor and Slot 


Design Details of Rotor and its Windings 

Length of Air gap(mm) 

lg 

i* 

0.8000 

Diameter of Rotor (mm) 

D 

D 

211.8 

Rotor Slot Pitch (mm) 

spr 

yss r 

17.0580 


Number of rotor slots that should not be used has been listed such that the slot combinations are suitable for stator 
and rotor for the given poles of the motor. The suitability is to avoid cogging, crawling and magnetic noise for the squirrel 
cage motor. Similar to air gap value the number of suitable rotor slots also may be optimised to get the best performance of 
the motor. This has not been done in the program. 


Table 26: Selection of Number of Rotor Slots 


Slots should not be equal 
to: 

36,44,40, 28 

Should not be equal to 

47,46,41,40 

No of Rotor Slots Selected 

39 
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XVIII STEP 3 OF PART III - ROTOR COPPER LOSSES 

In this step equivalent rotor current. Size of rotor bars are calculated. 

• Equivalent Rotor Current (I ph2 ) 

It is assumed that 85% of ampere turns get transferred to the rotor Ampere turns on stator = 3 I ph i N p 

S 2 So 

Ampere turns on rotor = I 2 bar — SoI 2bar — = 0.85 x 3 I ph i N ph i Or 

_ 0-85 x 6 l-phl Nph .1 


End ring current l 2endring = ^ X ^ 2 p bar 


• Size of Rotor Bars or Rotor Conductors 

Cross sectional area of rotor bar 


r* ^2 bar 

r c2bar _• 

Jbar 


Cross sectional area of rotor end ring 


r c2endring ~ ; V 

J endring 

• Design Calculations of Rotor Bars 

Selection of cross section of bars and end rings is a compromise between a) High rotor resistance which gives 
good starting torque and b) Low rotor resistance which gives low copper losses. Rotor bar (width to depth) geometry 
depends on the torque-speed characteristics and starting torque values. 

For a continuous service motor with high value of efficiency requirement low rotor resistance is aimed at. 
Rectangular shaped bars and slots are generally preferred in rotor. 

Current density in rotor bar is = 5 to 7 A/mm 2 for copper materials and for aluminium bar, the range is 2.2 to 4.5 
A/mm 2 and a higher limit value is chosen for motors analysis. 

Generally rotor bars and end rings are made of Aluminum material and stator conductors are made of copper 
material. It may be observed that the values are higher than stator conductors because rotor conductors are bare that is no 
insulation so better heat conduction resulting in better cooling so j is more. Current density in end ring is same as current 
density in bar. 

Table 27: Design Calculations of Rotor Bars and End Rings 
Step 3 of Part III 


Equivalent Rotor Current! A) 

Ir 

Iph2 

27.01 

Rotor bar Current! A) 

lb 

1 2 bar 

476.2 

Rotor bar CS(mm A 2) 

Ab 

Fc2bar 

82.3200 

Bar (w*t)mm 

a b (w*t) 

F c2bar 

(w*t) 

14. OX 6.0 

Rotor Slot 

sw.sh 

(w*h) 

14.5x6.5 

Length of Bar (m) 

Lb 

L 

257.000 
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Resistance/bar(m.ohm) 

Rb 

R b 

0.0656 

Losses in Rotor Bars (W) 

Pcub 

D 

1 cub 

579.8 

End ring Current 

Ie 

^2 endrin. 

1477.9 

Area of end ring (mm A 2) 

Ae 

Fc2endri 

246.3 

Resistance of end ring(m.ohm) 

Re 

Re 

0.0433 


Select current density in rotor bars and end rings and from rotor bar and end ring currents get their cross sectional 

areas. 


Calculate rotor bar and end ring resistances and hence copper losses. Equivalent rotor resistance is equal to total 
rotor copper losses divided by rotor current. Value of slip at full load is determined by rotor resistance. 

Resistance of one bar = 0.021 x 10" 6 — - — 

Fc2bar 

Cu loss in bars = S 2 x /f bar x Resistance of end ring 
- 0 0° 1 X lQ- 6 ^ n( ^ D ~ 2 9- 2 dibar ) 

r c2 end ring 

Cu loss in end rings = 2 x l£ en d r ing x Resistance of one end ring 
Total cu loss = Cu loss in bars + Cu loss in end rings 

Table 28: Rotor Copper Losses - Step 3 of Part III 


Used for 

MATLAB 

variable 

Variable 
used in 
Equations 

Value 

Bars 

Pcub 

Pcub 

579.8 

End. rings 

Pendring 

Pendring 

189.3 

Total 

Pcu,r 

P 

x cu,r 

769.1 

Equivalent 
Rotor res 'Q 

Rr 

R r or R 2 

0.351 


(23) 

(24) 


XXIII STEP 4 OF PART III - PARAMETERS AT 1/3 POSITION FROM NARROW END 

• Flux Density in Rotor Tooth and Core 

Note: The value of flux density in rotor tooth and core is same as flux density in stator tooth 

-j rd 

Diameter of rotor at — of tooth height from narrow end 
Di = D -2 lg - - h t2 x 2 

t-h t2 b 3 12 

End ring depth = Rotor bar (slot depth) + clearance 

-jTd 

Slot pitch at — of tooth height from narrow end 

D l 
t ? h t 

x n 

sg2-h t2 s 2 

From Narrow End 


(25) 


(26) 


2 rd 

Width of the tooth at — of tooth height 
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-1, — T _i, — bc 2 

t2 ~h t2 Sg2-h t2 sz 


Area of one stator tooth at — of tooth height = b „ u x K ) l 

3 b t2-h t2 1 


Area of all the stator teeth under one pole 


4 i, = Area of one tooth x No of teeth per pole 

t2.-h. t 2 


= b 1 xKjlx — 

t2-h t2 1 P 


n{o-2S~h t2 x2\ 


x K t lx — 


So mean flux density in teeth 


<Pi 


B t2\h t2 A t 


t^h t 2 


(27) 

(28) 


(29) 


(30) 


■^rcl 

Details of the rotor at — depth from the narrow end of rotor tooth section 


^rd 

Length of magnetic path =Pitch of rotor at — depth from minimum dimension of bar = L m 
Tooth width at bottom T b 

B s = Rotor Tooth Average Flux density at 1/3 position 

XXIV PARAMETERS AT THE POSITION OF 1/3 HEIGHT FROM NARROW SLOT END 

Parameters at 1/3 height from narrow slot end are evaluated to evaluate the average flux density 


Table 29: Parameters at 1/3 Position from Narrow End of Slot 


Used for 

MATLAB 

Variable 

Variable 
used in 
Equations 

Value 

Rotor Slot-Pitch(mm) 
Equation (28) 

SPrl3 

T 1 
sg2^h t2 

17.0580 

1/3 Position 

Slot-Pitch(mm) 



15.50 

tooth width) mm) 

equation (29) 

Wtrl3 

b i 
t2^h t2 

9.00 

T ooth-Flux-Dens 
(1/3) (T) 

Brt 

B i 

t2 gfl£2 

0.9369 

Tooth-Flux-Dens- 

Max(T) 

Brtmax 

Brtmax 

1.4054 

(Permissible: 1.2 to 1.5) 

Depth of Rotor core 
(mm) 

dcr 

dc.r 

30.45 


XXV STEP 1 OF PART IV NO LOAD PARAMETERS 

Estimation of No load parameters especially the no load current and no load torque and ensuring that they are not 
very high is very important for the good design of motor. Aluminium is lighter and has good forming characteristics so that 
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aluminium conductors may be ensured to have perfect contact between the slots and conductors. Copper offers low 
electrical and thermal resistances but perfect thermal contact may not be ensured if the conductors are not to be moulded. 
Moulding temperature associated problems have to be handled to consider copper conductors for moulding. 

Table 30: Variables and their Equivalents of MATLAB -Part IV 


Part IV 

MATLAB 

Variable 

Used for 

Variable 
used in 
Equations 

PnL 

No load loss 

P„L 

Iw 

No load phase current 

Iw 

Stator 

ATSC 

Amp-turns for stator core 

ATSC 

PfO 

No load power factor 

PfO 

ATST 

Amp turn for Stator teeth 

At2htl 

Spitch 

Stator Slot Pitch 

Spitch 

atst 

Amp-turns/m for stator teeth 

at 2 htl 

ATS 

Total AT for stator 

ATS 

Rotor 

WssO 

Stator slot opening 

WssO 

Spitch 

Rotor slot pitch 

Spitch 

kOl 

Carter coefficient 

koi 

K02 

Carter coefficient 

K„2 

SprO 

Rotor slot pitch near air gap 

SprO 

WsrO 

Rotor slot opening 

Wo, 

kgr 

Gap coefficient for Rotor 
Slots 

Kc2 

kgs 

Gap coefficient for Stator 
Slots 

Kci 

kg 

Air gap coefficient 

kg 

Lgd 

Effective air gap 

V 

Aag 

Air gap area/pole 

A ag 

Bg 

Flux density in the air gap 
(Bg) 

Bg (37) 

B30d 

Gap flux density at 30 from 
the centre of pole 

B 30 

Btr30 

Gap flux density in rotor 
tooth at 30 from the centre of 
pole 

Btr30 (38) 

atrt 

Amp-turns/m for Rotor tooth 

atr.t 

ATRT 

Amp-turns for Rotor tooth 

At2ht2 

atrc 

Amp-turns/m for Rotor core 

ak.v 

ATRC 

Amp-turns for Rotor core 

AT, V 

Im 

Magnetising current 

I m 

10 

No load phase current 

10 


The performance of the motor without load or starting of the motor is given by the following table. 

Table 31: No Load Parameters 


Used for 

MATLAB 

variable 

Variable 
used in 
Equations 

Value 

NO-Load Losses (W) 

PnL 

P„i 

1636.7 


www.tjprc.org 


editor@tjprc. org 



96 


A. Ravi Prasad & K Prahlada Rao 


Table 31: Contd 


Wattful Current (A) 
or core loss current 

Iw 

Iw or I c 

1.369 


Various values of Limits in the design of motor are listed in the following table. 


Table 32: Limiting Flux Density Values of Motor Components 


Position 

Typical Flux Density 
Range (Ref. -Say 
Author) 

Maximum Flux 
Density (Ref -Lipo 
Author) 

Air gap B s 

0.65 -0.82 T (average) 


Stator yoke 

1.1 -1.45 T (peak) 

1.7 T 

Stator teeth 

1.4 -1.7 T 

2.1 T 

Rotor yoke 

1.2 T 

1.7 T 

Rotor teeth 

1.5-1. 8 T 

2.2 T 


XXVI. ELECTRICAL AND MAGNETIC PARAMETERS 


The effect of change in magnetic flux due to the presence of slots on the laminations of stator and rotor empirical 
factors K cl> K c2 ,T sal =, i sgl = are eva | uatec | anc j 


Where K cl = Gap contraction factor for stator 


K cl 


L sg i 


L sg l 


■ (31) & r sgl = 


n D 
Si 


(32) 


The values of ratio of stator slot opening and air gap for stator and rotor are found out. Using the values of Slot 
pitch for stator and rotor the values of gap contraction for stator (32) and rotor (34) are calculated. Carter’s coefficients 
both for Stator and Rotor are found out from the interpolations. 


Carter coefficients for stator and rotor are obtained from the ratios: a) Stator slot width to air gap length ratio and 
b) Rotor slot width to air gap length ratio 


Kc 2 : 


l sg 2 


L sg2~ 


■ (33) & T sg2 = 


n(D-2lg ) 


(34) 


and Kqi & K 02 are constants or Carter coefficients for rotor and stator respectively 


IC, and k c2 are gap contraction factor for stator and rotor respectively 


l g ' = overall gap contraction factor = K c i x k c2 (35) 

Effective air gap length l g ' = l g K c iK c2 (36) 


Figure 1 



XXVII. MMF REQUIRED IN AIR GAP AT S 

H = — B AT/m (37) 

Mo 
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H = ~ #3o°l g ’ (38) 

Mo 

H = — B 30 o l g K cl K c2 (39) 

Mo 

XXVIII. FLUX DENSITY DISTRIBUTION 

Flux density at 30° from direct axis 
= flux density at 60 H from inter-polar axis 
So B° 0 = B° 0 = B ml Cos 30° 

= |B”y=1.36B” (40) 

For all practical purposes this value is modified to 

B 3 ° 0 =1.35B” (41) 

XXIX. STEP 3 OF PART IV 

For the calculation of Ampere -turns for Air gap. Stator tooth. Stator core. Rotor tooth, and Total No-load AT, 
Magnetizing current. No -Load PF. 


Table 33: Carters’ and AT/m Coefficients 


Interpolated Values Of Caters’ Coefficients 

Carter’s coeff - 
Stator 

kOl 

koi 

0.600 

Carter’s coeff - 
Rotor 

k02 

k()2 

0.395 

Stator Core 
AT/m 



446.7 

Teeth AT/m 



986.0 


XXX. STEP 4 AND STEP 5 OF PART IV 

In these two steps, the flux densities at locations of stator yoke, stator tooth, rotor tooth and rotor core (yoke) and 
air gap are calculate first. Then corresponding Ampere turns/ Unit length of magnetic flux path (AT/m) are calculated. 
With flux path length known, the values of AT for all locations and their sum are calculated described in steps of a) to f). 

Table 34: Length of Path and Ampere Turns at Locations -Formulas 


S .No 

Part 

Length 
of Path 

Flux 

Density 

at 

(AT/m) 

ATp 0 i 

e-Pair 

1 

Stator 

Yoke 

ly 

By 

at y 

ATy 

2 

Stator 

Tooth 

2h tl 

B ', 

t — h.i 

3 

a f2htl 

AT 2 hti 

3 

Air Gap 

2L’ 

b 30 o 

«lt2g’ 

AT 2r - 

4 

Rotor 

Tooth 

2h t2 

B 

3 rZ 

a Lht2 

AT 2ht2 

5 

Rotor 

Yoke 

'ry 

B ry 

atry 

AT ry 
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AT 


pole-pair 


AT X = £ 


Magnetising mmf is calculated with the help of magnetising curve and physical dimensions of stator teeth, stator 
core, rotor teeth, rotor core and air gap. 

• Length of path for air gap (21 g ’), flux density for air gap B 30 °, at 30(AT/m) -at 2< ; and AT 2o - for one pole are 
calculated 


i —n t j 

• Length of path for stator teeth (2h,i), flux density for stator teeth, 3 at 2 hti at (AT/m) and AT for stator teeth 
are calculated 

• Length of path for stator core or yoke (l y ), flux density for stator core (B y ), at (AT/m) at y and AT y for stator core a 
(AT 2 hti) are calculated 

• Length of path for rotor teeth (2h t2 ), flux density for rotor teeth, at (AT/m) and AT for rotor teeth are 

• calculated 

• Length of path for rotor core or yoke (lrj,), flux density for stator core (B ry ), at (AT/m) at ry and AT ry for stator core 
a (AT ry ) are calculated 

Table 35: Length of Path and Ampere Turns at Locations -Values 


Used for 

MATLAB 

Variable used 

Value 

variable 

in Equations 

Core 

ATSC 


38.5 

Teeth 

ATST 


42.4 

TOTAL 

ATS 


80.9 

ROTOR AT 

Core 

ATRC 


17.8 

Teeth 

ATRT 


5.5 

Table XXXV: Contd., 

TOTAL 

ATR 


23.3 

TOTAL AT 

Stator 

ATS 


80.9 

Rotor 

ATR 


23.3 

Airgap 



496.0 

TOTAL 

ATT 


600.2 


XXXI. STEP 1 OF PART V - NO LOAD MOTOR PERFORMANCE 

Determination of Core Loss Current (I c ) 

Total iron losses P, = Iron loss in teeth + Iron loss in yoke 


Active/ Wattful component of no load /Phase current I c = — l - — 

F 3 v phl 


(42) 


The weight of iron in the stator is found and iron losses are calculated from material data sheets of losses in W/kg 
depending on flux density and frequency. Iron loss in rotor is negligible because of low frequency and this loss is generally 
accounted as part of stray load loss or additional loss. 
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Magnetising Current 


^30 


P 

2 


1 


l-17Nphl Kdl 


(43) 


Calculation of No Load Current (I 0 ) 

To solve the magnetic circuit of motor, the no load current has to be determined. Thus one has to estimate the 
magnetizing component (Im) and then core loss component (I c ) of current. 

No load current I 0 




No load power factor Costp 0 = p 


Table 36: Performance of Motor at No Load 


No-Load Current (A) 

Eq. 

Core loss 

Iw 

Iw 

1.37 

(42) 

Magnetising 

Im 

Im 

8.95 

(43) 

No Load 

Io 

Io 

9.06 

(44) 

Power factor 

Pf 

Pf 

0.151 

(45) 

Ratio of No Load Current to Phase Current 


Io/Iph 

ratio 

Io/Iph 

ratio 

0.285 



(44) 

(45) 


The six impedances are stator resistance Ri, stator leakage reactance X 0 i, magnetizing reactance X 0m , core loss 
Resistance R m , rotor leakage reactance X a2 , and rotor resistance R 2 . 

XXXII. STEP 3 OF PART V - SHORT CIRCUIT RESISTANCES -MOTOR PERFORMANCE 

MATLAB program is structured such that there is a scope for introducing a loop which can take care of different 
geometries of slots after the dimensional details of slot geometries are input for carrying out the optimisation. Generally 
this way of achieving motor optimisation is considered as the difficult strategy and optimisation is not done from this 
aspect. 


Specific Slot Permeability 

Specific performance is evaluated based on 

Specific permeability for a tapered stator slot is 


^sp,s — Bo [ 


2hi 


2h? 


h 3 h 4 -j 

■ 1 ] 

(3W 2 +W ss ) (Wi+W 2 ) (Wi+W 0 ) Wq 

Specific permeability for a tapered rotor slot is 


* r hi h 2 2 h 3 h 4 

^sp,r = BO [— + — + I- ] 

•3W w (w+Wq ) Wo 


(46) 


(47) 


www.tjprc.org 


editor@tjprc. org 


100 


A. Ravi Prasad & K Prahlada Rao 


Total Specific Slot Permeability 

X s = (Specific slot permeability for stator slot+ specific slot permeability for rotor slot). 


— Ap.s t-sp.r 

Slot Leakage Reactance 

8 jt f Nj?,, L 


Xs = ■ 


N ph 1 


P qs 


Overhang Leakage Reactance 


Xa _ 


0.64 tt f Np h n 0 T 
" P QsYss 


TT 7T D , TT D 

Here, r =— and y ss = — 


Zigzag Leakage Reactance 


Core reactance X m = 


^ph 

Im 


v _ 5 X m r 1 . l-i 

Xz -6 


Here, rotor slot per pole per phase, q, = 


Sr 

3p 


(48) 


(49) 


(50) 


(51) 

(52) 


Total Leakage Reactance per Phase 

X = X, = (X S +X 0 +X Z ). 

(53) 

z= Vr 2 + x 2 

(54) 

Short Circuit Current 

T V ph 

isc- z 

(55) 

Pf -- 

risc z 

(56) 

Overload capacity = — 

>ph 

(57) 


Table 37: Total Reactance (Short Circuit Resistance) 


Total Reactance Equation 

Slot 

(47) 

x s 

Xs 

0.925 

Overhang 

(48) 

X 0 

xo 

0.909 

Zig-Zag 

(49) 

x z 

Xz 

0.648 

SHORT 

CIRCUIT 

Resistance 


R 

R 

0.67 

Total leakage 

REACTANCE 

(X) 

(51) 

X 

X 

2.481 
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z 

(52) 

Z 

Z 

2.57 

Iso 

(53) 

Iso 

Isc 

155.0 

pf 

(54) 

pfsc 

pfsc 

0.260 

( Isc / Iph ) Of ( Isc 

/Ifl) 

(55) 

Iratio 

Iratio 

4.879 


XXXIII. STEP 4 OF PART V LOAD LOSSES 

Load losses are the losses at the operating load losses will finally evaluate the performance of the motor. For some 
motors which operate at different loads performance at different loads also are evaluated. 

Load loss= stator copper loss + rotor copper loss + No load 


Table 38: Losses of Motor at Operating Load 


No Load Losses 

PnL 

P„r 

961.9 

Stator copper losses 

Pcus 

Pcus 

769.1 

Rotor copper losses 

Pcur 

p 

A cur 

1636.7 

Total Losses 

Pt 

P total 

3367.7 


XXXIV. STEP 5 OF PART V - TOTAL WEIGHT AND KG/KW 

Estimation of weights of all the components and total weight of motor of a particular capacity of motor is essential 
where optimization leading to a lighter motor for the given capacity motor is a requirement. Total weight of all the main 
components and specific power is calculated. 


Table 39: Weights of Motor Components 


Total Weight Of The Components 

Stator copper 

wcus 

w cus 

27.4 

Stator teeth 

wt 

w t 

19.2 

Stator core 

wsi 

W si 

46.8 

Rotor core 

wri 

W ri 

32.1 

Rotor bars 

wcur 

w cur 

7.3 

End rings 

wcue 

w cue 

2.2 

Total 

wtotal 

Wtotal 

136.5 

Kg/KW 

kgpkw 

Wtota/KW 

4.55 


XXXIV. STEP 6 OF PART V -MOTOR PERFORMANCE LOAD 

In this step motor performance parameters namely efficiency. Slip at FL, Starting torque. Maximum Output (KW) 
and Temperature rise are estimated. 

q = output / (output+ No load loss + Load loss) 

R inp = Input power to Rotor = KW of motor + Rotor copper losses + Losses due to friction and windage 
Slip at FL = Rotor Copper losses / Input power to Rotor 


Table 40: Motor Performance Parameters 


Efficiency % 

EFF 

n 

89.91 

Slip at FL % 

SFL 

Sfl 

2.475 

Starting Torque 

Tst 

T st 

0.816 

Maximum Output 

Pmax 

P 

A max 

69.2 
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XXXV. STEP 7 OF PART VAPPROXIMATE TEMPERATURE RISE ESTIMATION 

Highly approximate method of estimation of temperature rise of stator frame has been mentioned in this step. In 
this evaluation end caps are considered as part of the frame and the temperature corresponds to broad average temperature 
rise of the frame which includes end caps on both sides of TEFC motor. 

For all the motors efficiency (r|) of fan is considered as 0.5 

Inner area of stator bore = A inner =(ji x D x L*2.5 + (2.* n *(D+50)*0.04)/le06 

Peripheral velocity = v = 7ixDxnsx efficiency of fan 

Effective area factor = Af actor = (1 + 0.1 v) 

Effective area = A efanner = A factor x A inner 

Outer stator area = A outer = n x D out x L 

Total dissipating area = Acool_total = A eff j nne . r + A outer 

Tr = 0.03 *Pst/Acool_total; 


Table 41: Approximate Temperature Rise 


Inner 
area of 
Stator 
Bore m 2 

Peripheral 

Velocity 

m/sec 

Effective 

Area 

Factor 

m 2 

Effective 
Area m 2 

Outer 

Stator 

Area 

m 2 

Total 

Area 

m 2 

Temperature 
Rise of 
Frame 
Surface 

0.346 

8.37 

1.837 

0.926 

0.234 

1.16 

59.45 


XXXVI OPTIMISATION OF MOTOR 


For optimization of motor current density values, slot height to width and slot fill factors are varied for round 
conductors and current density values, height to thickness of bar and height to width of slot are varied. The results are 
presented in the following tables. J and Thick refer to current density and insulation thickness in the following tables. 


Table 42: Optimization of Motor Using Rod Conductors 


Rod Conductors are Used 

WXH 

Rt 

DL max 

T r 

D 0 

h 

Kgp 

Kw 

IObyl 

Rods are 3.2 mm dia x 0.5 thick and J 

is 3.5 A/mm 2 

7.0X21 

1.31 

45.7 

320 

91.67 

4.02 

0.30 

6.3X19 

1.21 

44.7 

320 

91.78 

3.99 

0.30 

6.1X24 

1.15 

46.3 

330 

91.58 

4.18 

0.30 

5.4X22 

1.09 

45.2 

320 

91.73 

4.12 

0.30 

Rods are 3.0 mm dia x 0.5 thick and J 

is 3.5 A/mm 2 

6.8X20 

1.28 

48.4 

320 

91.39 

3.90 

0.30 

6.1X18 

1.19 

47.4 

320 

91.49 

3.87 

0.30 

5.9X24 

1.13 

48.9 

330 

91.30 

4.06 

0.30 

5.3X21 

1.07 

47.9 

320 

91.44 

3.99 

0.30 

Rods are 2.8 mm dia x 0.5 thick and J is 3.5 A/mm 2 

6.7X20 

1.26 

51.1 

320 

91.11 

3.81 

0.30 

6.0X18 

1.18 

50.1 

310 

91.24 

3.76 

0.30 

5.8X23 

1.12 

51.6 

330 

91.02 

3.96 

0.30 

5.2X21 

1.06 

50.7 

320 

91.15 

3.90 

0.30 
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%M 

ax. 

Eff. 

Min. Kg 
/KW 

Min Temp- 
Rise 

Min 

IOIr 

ratio 


91.78 

3.76 

44.70 

0.30 

Variant 

2 

0 

2 

10 


Table 43: Optimization of Motor Using Bar Conductors 


Bar conductors are used 

WXH 

® tmax 

T r 

Do 

il 

Kgp 

Kw 

IObyl 

Bar C/s is 12 x 6.5 Current Density 3.5 A/mm2 

8.7X48 

1.39 

55.6 

370 

90.47 

4.70 

0.31 

8.7X48 

1.39 

55.6 

370 

90.47 

4.70 

0.31 

Bar C/s is 1.8 x 5.0 and Current Density is 3.5 A/mm2 

10.2X40 

1.85 

53.4 

360 

90.74 

4.40 

0.39 

10.2X40 

1.85 

53.4 

360 

90.74 

4.40 

0.39 

Bar C/s 

is 1.2 x (i 

.0 and Current Density 

is 4.0 A/mm2 

8.7X46 

1.41 

56.6 

370 

90.37 

4.56 

0.31 

8.7X46 

1.41 

56.6 

370 

90.37 

4.56 

0.31 

Bar C/s is 1.8x 4 

.5 and Current Density is 4.0 A/mm2 

10.2X38 

1.89 

54.7 

350 

90.63 

4.22 

0.41 

10.2X38 

1.89 

54.7 

350 

90.63 

4.22 

0.41 

Bar C/s is 1.2x 5.5 and Current Density is 4.5 A/mm2 

8.7X43 

1.43 

57.9 

360 

90.27 

4.41 

0.31 

8.7X43 

1.43 

57.9 

360 

90.27 

4.41 

0.31 

Bar C/s is 1.8x 4 

.0 and Current Density is 4.5 A/mm2 

10.2x36 

1.93 

56.6 

350 

90.44 

4.07 

0.42 

10.2x36 

1.93 

56.6 

350 

90.44 

4.07 

0.42 


Max. 
Eff. % 

Min. Kg 
/KW 

Min A0 

Min IOIr 
ratio 


90.74 

4.07 

53.42 

0.31 

Variant 

4 

0 

4 

10 
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APPENDICES 


%Annexure 

MATLAB Program for Traditional Design and 
% Optimization of Performance of 30 KW SCIM Motor 
f2=fopen('motor_30KW_m','w'); 

%<— Step I of part I —Basic Dimensions of Motor-> 
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optim=l ;cond_type=l ;slot_fill_factor=0.4; 

D=213; L=207;KW=30.0;f=50;P=4;%Specification/lnput Data > 

Ns=120*f/P;ns=Ns/60;Kw=0.955;insH=6;insS=0.5; %Initial Values 
KWinp=30; 

% Step 2 of Part I - Electro magnetic parameters of motor 
V=398.37;ki=0.96; Tstrip=1.50 ; 

insW=3.4;Hw=4;HL=l;insH=6;ki=0.92;Bc=1.35;Vph=V;%Initial Values 
Zr=l;kwr=l;cdb=6;Tb=6;cde=6;dd=0.05;Brc=1.35; 

sn=0;M 1 =0;M2=0;M3=0;M4=0;EFFmax=90;minKgPkw=9 ;minTr=75 ;minIObyI=5 ; 
Vph=V; P=4; spp=4;Zsw=3; 

% Step 3 and step 4 are for checking the limits like bearings and frames 
%Initial Values -Standard Curves/Tables for%Data — > 

SD=[0.1 0.15 0.2 0.3 0.4 0.5 0.75 1]; 

SCDSW=[4 3.8 3.6 3.5 3.5 3.5 3.5 3.5]; 

B=[0.8 1.2 1.6 2 2.4]; WpKg=[7 15 24 34 50]; 

%Step 6 of Part I of traditional design 
SKW=[1 2 5 10 20 50 100 500]; 

SBav=[0.35 0.38 0.42 0.46 0.48 0.50 0.51 0.53]; 

Sq=[16e3 19e3 23e3 25e3 26e3 29e3 31e3 33e3]; 

SKWa=[5 10 20 50 100 200 500]; 

SPF4P=[0.85 0.86 0.88 0.9 0.91 .92 ,93];%Table: for 1500RPM 
SEFF4P=[.85 .87 .88 .9 .91 .93 ,94];%Table for 1500RPM 
SD=[0.1 0.15 0.2 0.3 0.4 0.5 0.75 1]; 

SCDSW=[4 3.8 3.6 3.5 3.5 3.5 3.5 3.5]; 

% BH Curve for— 0.5mm, LOHYS Quality > 

BBone= [ .1 .2 .3 .4.5 .6 .7 .8 .9 1 ]; 

BBtwo= [1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2]; 

Hone= [50 65 70 80 90 100 110 120 150 180] ; 

Htwo= [220 295 400 580 1000 1400 5000 8900 15000 24000]; 
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BB= [BBone BBtwo]; 

H = [Hone Htwo]; 

% Carters Coefft for Air Gap— > 

Ratio=[0 1 2 3 4 5 6 7 8 9 10 11 12]; 

CC= [0 .18 .33 .45 .53 .6 .66 .71 .75 .79 .82 .86 .89] ; 

CC1= [0 .14 .27 .37 .44 .5 .54 .58 .62 .65 .68 .69 .7]; 

%Step 5 of Part I of traditional design 
pf=interp 1 (SKWa,SPF4P,KW, 'spline'); 
eff_data=interp 1 (SKWa,SEFF4P,KW, 'spline'); 

Bav=interpl(SKW,SBav,KW, 'spline'); 
q=interpl(SKW,Sq,KW, ’spline’); 

Iph=KWinp* 1 e3/(3 * Vph*pf); 

CO=ll*Kw*Bav*q*eff_data*pf*le-3; % Equation (4); 

KWinp=KW/eff_data; 

Fs=F;Fi= ki*Fs; PP=pi*D/P; FbyPP=F/PP 
% if FbyPP <0.8 II FbyPP >2 continue; end; 
v=pi*D*ns/1000 
if v >30 continue;end; 

FI=pi*D/P*F*Bav/le6 ; 

Tphi=Vph/(4.44*f*FI*Kw);S=spp*P*3; 

CDSW_int=3.5; CDSW_end=4.5;CDSW_ave=(CDSW_int+ CDSW_end)/2; 
if(cond_type==l) % Circular conductors 

F2_int=3.0;F2_end=4.0;F2_incre=l;% loop two is variation of slot height and slot width ratio 

F3_int=0.4;F3_end=0.5;F3_incre=0.1;% loop three is varying slot_fill factor 

end; 

if(cond_type==2)% Bar conductors 
F2_int=1.25; F2_end=1.75; F2_incre=0.5;% Bar thickness 
F3_int=2;L3_end=3;F3_incre=l; % Bar height to width ratio 
T strip=( F2_int+F2_end)/2 ; 
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end; 

L2_ave=(L2_int+L2_end)/2; 

L3_ave=(L3_int+L3_end)/2; 

if optim==0 % Fixed values of current density and slot dimensions 
CDSW_int=CDSW_end;L2_int=L2_end; 

L3_int=L3_end; 

end; 

for CDSW =CDSW_int:0.5:CDSW_end; 

Asl=Iph/CDSW; 

if optim==0 Asl=Iph/CDSW_ave;end; 

SPitch=pi*D/S; 

% if SPitch <18 II SPitch>=25 continue; end; 

Zphi=2*Tphi; sph=S/3; Zsl=Zphi/sph;Zs=ceil (Zsl); Tph=Zs*sph/2; 

slots_per_phase=sph;cond_per_phase=Zphi; 

cond_per_slot= (cond_per_phase/slots_per_phase); 

FI=Vph/(4.44*f*Tph*Kw); 

if optim==0 L2_incre=l;L3_incre=l;end; 

for L2_value =L2_int:L2_incre:L2_end; 

if(cond_type==2) 

Tstrip=L2_ave; 

if optim==l Tstrip=L2_value;end; 

Hstrip 1 = As 1 /L2_value; 

FIstrip=ceil(Hstripl*2)/2 ; 

WbyT =Hstrip/T strip 

% if WbyT <2.5 II WbyT >=3.5 continue; end; 

As=0.967*Hstrip*Tstrip; 

end; 

for L3_value=L3_int : L3_incre: L3_end ; 
if(cond_type==2) 
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zsw=L3_value; 

Ws=(Zsw*(Tstrip+insS)+insW); 

Zsh=Zs/Zsw; 

Hs=(Zsh*(Hstrip+insS)+Hw+HL+insH+2);% Hw is height of wedge and 2 is clearance 
end; 

if(cond_type== 1 ) 

As=Iph/CDSW; 
cond_dia = (As*4.0/pi) A 0.5; 

Htby W d=L2_ave ; 
slot_fill_factor=L3_ave ; 

if(optim==l) HtbyWd=L2_value;slot_fill_factor=L3_value;end; 
slot_area= (cond_dia+2*insS) *cond_per_slot/slot_fill_factor; 

Slot_wide =(slot_area /HtbyWd) A 0.5;Slot_hight= Slot_wide*HtbyWd; 

cond_along_width=(Slot_wide/(cond_dia+2*insS)); 

cond_along_hight=(Slot_hight/(cond_dia+2*insS)); 

Ws=Slot_wide;Hs=Slot_hight; 

Zsh=cond_along_hight;Zsw=cond_along_width; 

Hstrip=cond_dia+2.0*insS; 

end; 

%Step 3 of Part II 

D13=D+2/3*Hs; 

spl3=pi*D13/S; 

Wtl3=spl3-Ws; 

B13=FI*P*le6/(Li*Wtl3*S);% Equation 15 
Btmax=1.5*B13; 

Lmt=(2 *L+2 . 3 *PP+240)/ 1 000; % Equation 16 

Rph=0.021*Lmt*Tph/As;% Equation 17 

Pcus=3 *Iph A 2*Rph; Wcus=Lmt*Tph*3*As*8.9e-3; Flc=FI/2; 

Ac=Flc * 1 e6/Bc ;Hc= Ac/Li ; D01=D+2*(Hs+Hc); DO=ceil(D01/10)*10; 
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Hc=(DO-D)/2-Hs; 

% Table XXIII Iron losses and weight of stator 
PitpKg=interpl(B,WpKg,Btmax, 'spline'); 

PicpKg=interpl(B,WpKg,Bc, ’spline'); 

Wt=Li*Wtl3*S*Hs*7.8e-6; %Weight of teeth 
Dmcs=D+2*Hs+Hc; Wc=Ac*pi*Dmcs*7.8e-6; %Weight of core 
Pit=PitpKg*Wt;Pic=PicpKg*Wc; 

% ROTOR > 

Zr=l;kwr=l;cdb=6;;cde=6;dd=50;Brc=1.35; % Initial Values 
kws=Kw;Ss=S;Lgl=0.2+2*sqrt(D*L/le6); 

Lg=ceil(Lg 1*1 00)/ 1 00; 

% Rotor slot numbers that need to be avoided for preventing rotor locking 
dl=Ss-3*P;d2=Ss-P;d3=Ss-2*P; d4=Ss-5*P;d5=Ss-l;d6=Ss-2;d7=Ss-7;d8=Ss-8; 
n_lams= 1 ; 

LG= [0.8 0.8 0.8 0.8];WSS0=[4.0 4.0 4.0 4.0];WSR0=[2.0 2.0 2.0 2.0]; 

H2=[1.6 1.6 1.6 1.6];H2R=[0 0 0 0]; H3R=[0 0 0 0]; H4R= [0.5 0.5 0.5 0.5]; 
for L4 =l:n_lams 
Lg=LG(L4);% Lg=LG(L4); 

Dr=D-2 *Lg; Sr=39 ; Tb=6.0; 

%Step 1 and Step 2 of Part III 
sp2=pi*Dr/Sr; Ir=0.85*Iph; 

%Step 3 of Part III 

Ib=Ir*kws * S s*Zs/( kwr* Sr*Zr) ; %Equation 1 9 
Abi=Ib/cdb; Wb=ceil(Abi/Tb); 

Ab=Tb*Wb*0.98;Wsr=Tb+0.5;Hsr=Wb+0.5;Lb=L+50;Rb=0.021*Lb/le3/Ab; 
Dme=Dr-dd; Pcub=Ib A 2*Rb*Sr; Ie=Ib*Sr/P/pi; Ae=Ie/cde; 
Lme=pi*Dme/1000;Re=0.021*Lme/Ae; 
Pcue=2*Ie A 2*Re;Pcur=Pcub+Pcue;Rr=Pcur/(3*Ir A 2); 
Drl3=Dr-2*2/3*Hsr;sprl3=pi*Drl3/Sr 
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Wtrl3=sprl3-Wsr;Atr=Wtrl3*Li*Sr/P;Brt=FI*le6/Atr; 
Brtmax=Brt*1.5;Ac=FI*le6/2/Brc;dcr=Ac/Li; 
Pfw=0.01*KW*le3;PnL=Pit+Pic+Pfw;Iw=PnL/3/V; 
Wcur=Lb*Sr*Ab*8.9e-6;Wcue=Lme*2*Ae*8.9e-3; 
%(4)<— AmpTurns and Magnetizing-Current— Part IV— > 
W sso=WSS0(L4);W sro=W SR0(L4); %Initial Values 
atsc=interpl(BB,H,Bc, 'spline'); 

%Step 4 and Step 5 of Part IV 
Dcav=D+2*Hs+Hc ; ATSC=pi*Dcav/P/3 *atsc/ 1 e3 ; 
Bt30=B13*1.36;% Equation (41) for stator 
atst=interp 1 (BB ,H,Bt30, 'spline'); 

ATST=atst*Hs/1000; 

ATS=ATSC+ATST %Amp turns for tooth and core 
ratl=Wsso/Lg 

k01=interp 1 (Ratio, CC, rati , 'spline'); 
kgs=SPitch/(SPitch-Wsso*k01); %Equation 32 
rat2=Wsro/Lg; 

k02=interpl(Ratio,CC,rat2, 'spline'); %Equation 48 
spro=pi*Dr/Sr; %Equation 35 

kgr=spro/(spro-Wsro*k02); kg=kgs*kgr; Lgd=Lg*kg; 

% rat3=bvd/Lg; kv=interpl(Ratio,CCl,rat3, 'spline'); 

% if rat3 >=12 kv=0.7;end; 
Ld=L;Aag=pi*D/P*Ld;Bg=FI*le6/Aag; 

B30d=1.36*Bg; 

ATg= 0.796*B30d*Lgd*le3 
Btr30=Brt*1.36 % Equation 41 for rotor 
atrt=interpl(BB,H,Btr30, 'spline') % at of rotor teeth 
ATRT=atrt*Hsr/le3; 

Dcrav=Dr-2*Hsr-dcr; 
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atrc=interpl(BB,H,Brc, 'spline') 

ATRC=pi*Dcrav/le3/P/3*atrc; 

ATR=ATRC+ATRT % AT of rotor core and rotor teeth 

ATT=ATS+ATR+ATg; % Step 1 of part V 

Im=P/2*ATT/(1.17*Kw*Tph); %Equation 43 

I0=sqrt(Iw A 2+Im A 2);%Equation 44 

PfO=Iw/IO; I0byl=10/Iph; %Equation 45 

% (5) <— Step 3 of part 5— Short-Circuit-Current— >'); 

h2=H2(L4) ;h2r=H2R(L4) ;h3r=H3R(L4) ;h4r=H4R(L4) ; ks=l; %Initial Values 

if(cond_type== 1 ) Hstrip=cond_dia+2*insS;end; 

hl=Zsh*(Hstrip+insS);h3=Hw;h4=HL;bs=Ws; 

bO=Wsso;Lmdss=hl/3/bs +h2/bs +2*h3/(bs+bO) +h4/bO; 

hlr=Wb; br=Wsr; %Equation 45 

brO=Wsro;Lmdsr=hlr/3/br+h2r/br+2*h3r/(br+brO)+h4r/brO; 
Lmddsr=Kw A 2*S/Sr*Lmdsr; ssp=Lmdss+Lmddsr; gd=S/P/3; 
p=P/2;Xs=15.8*f*L*ssp*Tph A 2/(p*gd)*le-9 %Equation 47 
LOLmdO=ks *PP A 2/pi/SPitch/l 000; 

XO=15.8*f*LOLmdO*Tph A 2/(p*gd)*le-6 %Equation 48 
gs=S/P;gr=Sr/P;Xm=Vph/Im %Equation 49 
Xz=5/6 *Xm*( 1 /gs A 2+ 1 /gr A 2) %Equation 50 
X=Xs+X0+Xz; R=Rph+Rr;%Equation 51 
Z=sqrt(R A 2+X A 2); Isc=Vph/Z; 
pfsc=R/Z; RAT=Isc/Iph;%Equations 52,53 

%(6)< Performance >’); 

Pt=PnL+Pcus+Pcur; %Table XXXVII 
eff=KW/(KW+Pt/l 000)* 1 00 
if optim==l & eff<88.8 continue;end; 

Rinp=KW * 1 000+Pfw+Pcur; 

SFL=Pcur/Rinp * 1 00; Tst=(Isc/Ir) A 2*SFL/100; 
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Pmax=3 * Vph*(Isc-I0)/2/( 1 +pfsc)* 1 e-3 ; 
Acooll=(pi*D*(L*2.5)+2*pi*(D+50)*0.04)/le6; 

Acoo 1 2= Acool 1 *( 1 +0. 1 *v); 

Acool3=pi*DO*L/le6; 

AcoolT = Acoo 1 2+ Acoo 1 3 ; 

Pst=Pcus+Pit+Pic; 

Tr=0.03*Pst/AcoolT; 

% if optim==l & Tr>56 continue;end; 

Ars=W sr*Hsr*Sr;Dri=Dcrav-dcr; 

Wri=(pi*(Dr A 2-Dri A 2)/4-Ars)*L*7.8e-6; 

Wtot=1.01*(Wcus+W t+ W c+ Wri+ W cur+ W cue) ;KgPk w= W tot/KW ; 
if optim==l & KgPkw>6.1 continue;end; 

% End of TRaditional design Program — > 

if optim==0 

fprintf(f2, ’Design of 30KW,440V,50HZ,3-Ph SCIM\n'); 
fprintf(f2, ’\nlnput Data:’); 

fprintf(f2, ’\n '); 

fprintf(f2, ’\nParamter VALUES'); 

fprintf(f2, '\n— '); 

fprintf(f2, ’\nRating(KW) %5.1f,KW); 
fprintf(f2, ’\nVolts%5.0f,V); 
fprintf(f2, ’\nPoles%5.0f,P); 
fprintf(f2, ’\nHz %5.0f,f); 

fprintf(f2, '\nInterpolated values from curves: B a v = % 5 . 3 f , q = % 5 . 0 f , B a v , q ) ; 
fprintf(f2, ’eff=%7.2f,pl=%5.3f,eff,pf); 
fprintf(f2, '\nOutput Results:'); 
fprintf(f2, '\n==============='); 

fprintf(f2, ’\nParamter VALUES'); 
fprintf(f2, ’\n '); 
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fprintf(f2, ’\nOutput Coefft(CO)%6.2f ,CO); 

fprintf(f2, ’\nSync.Speed(rps) %5.2f,ns) ; 

fprintf(f2, '\nGross Length(mm) %5.1f,L); 

fprintf(f2, '\nNet iron Length(mm) % 5. If, Li) ; 

fprintf(f2, '\nStator Inner Dia(mm) %5.1f,D); 

fprintf(f2, 'VnPeriphoral Speed(m/s) %5.2f (Max. Limit: 30)', v) ; 

fprintf(f2, ’\nPole-Pitch(mm)%5.1f,PP); 

fprintf(f2, '\nL to PP ratio %6.4f (Around 1 >good)',LbyPP) ; 

fprintf(f2, ’\nSlots %5.0f,S); 

fprintf(f2, '\nSlot-Pitch( mm)%6.3f( Limit: 18to25mm)',SPitch); 
fprintf(f2, ’\nCond/Slot %5.0f ,Zs); 
fprintf(f2, ’\nTurns/Ph%5.0f,Tph); 
fprintf(f2, ’\nFlux/Pole(Wb)%6.5f,FI); 
fprintf(f2, '\nPhase Current (A) %5.2f,Iph); 

if cond_type==2 fprintf(f2, '\nBare Strip (w*t)mm %5.2fX%4.2f,Hstrip,Tstrip); 

fprintf(f2,’\nW and T Ratio %5.2f(Limit:2.5to3.5)',WbyT); 

end; 

if cond_type==l fprintf(f2, ’\nBare dia and insulation mm %5.2fX%4.2f,cond_dia,insS); 

fprintf(f2,'\nW and T Ratio %5.2f(Limit:2.5to3.5)',WbyT); 

end; 

fprintf(f2, ’\nArea of CS cond(mm A 2) %7.3f ,As); 

fprintf(f2, ’\nCt.density(A/mm A 2) %7.3f,CDSW); 

fprintf(f2, '\n values of W/KgTeeth=%5.2f Core=%5.2f,PitpKg, PicpKg); 

fprintf(f2, ’\n Strips (width*depth-wise) %4.0f X%2.0f ,Zsw,Zsh); 

fprintf(f2, ’\nSlot- Width (mm) %6.1f,Ws); 

fprintf(f2, ’\nSlot-Height(mm)Total %6.1f,Hs) ; 

fprintf(f2, ’\n d(l/3),SP(l/3),Wt(l/3)(m)%6.4f,%6.4f,%6.4f ,D13,spl3,Wtl3); 
fprintf(f2, ’\nSt-Tooth-Flux-Dens(l/3) %6.4f,B13); 

fprintf(f2, ’\nSt-Tooth-Flux-Density -Max(T) %6.4f(Limit:L6tol.8)',B13*L5); 
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fprintf(f2, '\nLength of mean-turn (m)%6.3f ,Lmt); 
fprintf(f2, "VnResistance/Ph (ohm) % 6.4f,Rph) ; 
fprintf(f2, '\ndepth of St.Core(mm) %6.2f ,Hc); 
fprintf(f2, '\nOuter Dia of St.Core(mm) %6.1f ,DO); 
fprintf(f2, ’\nSt.Cu.Loss(W) %6.1f,Pcus); 

fprintf(f2, ’\nWt of St-Teeth+Core(Kg)= %6.2f+%6.2f=%6.2f ,Wt,Wc,Wt+Wc); 
fprintf(f2, ’\nIron Loss=Teeth+Core(W)= %5.1f+%5.1f= %5.1f,Pit,Pic,Pit+Pic); 

fprintf(f2, ’\n-Part III— ROTOR- Design >’); 

fprintf(f2, '\nLength of Air-Gap(mm)%6.4f ,Lg); 
fprintf(f2, '\n Dia of Rotor(mm) ) % 6. If , Dr); 

fprintf(f2, ’\nNo of Slots Should NE to:%2.0f,%2.0f,%2.0f, % 2.0f, %2.0f,dl,d2,d3,d4); 

fprintf(f2, ’\nShould NE to:%2.0f,%2.0f,%2.0f,%2.0f,%2.0f ,d5,d6,d7,d8); 

fprintf(f2, '\nNo of Rotor Slots Selected %3.0f,Sr); 

fprintf(f2, ’\nRotor Slot-Pitch(mm) %6.4f ,sp2); 

fprintf(f2, ’\nEquivalent Rotor Ct(A) %6.2f ,Ir); 

fprintf(f2, ’\nRotor bar Ct(A) %6.1f,Ib); 

fprintf(f2, ’\nRotor bar CS(mm A 2) %6.4f ,Ab); 

fprintf(f2, ’\nBar (w*t)mm %5.1fX%4.1f ,Wb,Tb); 

fprintf(f2, ’\nRotor Slot (w*h)mm %5. lfX%4. If ,Wsr,Hsr); 

fprintf(f2, ’\nLength of Bar (m) %5.3f,Lb); 

fprintf(f2, ’\nResistance/bar(m.ohm) %6.4f ,Rb*le3); 

fprintf(f2, ’\nLosses in Rot. Bars (W) %6.1f,Pcub); 

fprintf(f2, ’\nEndring Ct(A) %6.1f,Ie); 

fprintf(f2, ’\nArea of end ring (mm A 2) %6.1f ,Ae); 

fprintf(f2, ’VnResistance of end ring(m.ohm) %6.4f,Re*le3); 

fprintf(f2, ’\nRot-Cu- Loss=Bars+End.rings(W)=%5. 1 f+%5. 1 1=%5. 1 f.Pcub.Pcue. Pcur); 
fprintf(f2, ’\nEquivalent Rotor res (Ohm) %6.3f ,Rr); 
fprintf(f2, ’\nRotor(l/3)Slot-Pitch(mm) %6.2f ,sprl3); 
fprintf(f2, ’\nRotor(l/3)tooth width(mm) %6.2f ,Wtrl3); 
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fprintf (f2, VRt-Tooth-Flux-Dens (1/3) (T) % 6.4f, Brt) ; 

fprintf(f2, ’\nRt-Tooth-Flux-Dens-Max(T) %6.4f(Limit:1.2tol.5)' ,Brt*1.5); 

fprintf(f2, '\nDepth of Rotor core (mm) % 5.2f ,dcr) ; 

fprintf) f2, '\n No-Load-Losses >'); 

fprintf) f2, ’\nNO-Load Losses (W) %5.1f,PnL); 
fprintf) f2, '\nNO-Load Wattful Current (A) %5.3f,Iw); 
fprintf) f2, ’\n Magnetizing-Current >'); 

fprintf(f2, '\nlnterpolated values of at/m of St:Core=%5.1fandTeeth=%5.1f,atsc,atst); 

fprintf(f2, '\nlnterpolated values of Caters Coeffts:k01=%5.3f,k02=%5.3f,kv= %5.3f, k01,k02); 

fprintf(f2, ’\nStator %5.1f+%5.1f=%5.1f ,ATSC,ATST,ATS); 

fprintf) f2, ’\nRotor AT :Core+Teetn: %5.1f+%5.1f=%5.1f .ATRC.ATRT.ATR); 

fprintf(f2, ’\nTotal AT :Stator+Rotor+Airgap: %5.1f+%5.1f+%5.1f= %5.1f .ATS.ATR, ATg.ATT); 

fprintf(f2, '\nNO-Load Current (A) Iw=%5.2f,Im=%5.2f and IO=%5.2f at pf=%5.3f,Iw,Im,I0,Pf0); 

fprintf(f2, ’\nIO/lph ratio %5.3f,I0byI); 

fprintf) f2, '\n Short-Circuit-Current >'); 

fprintf(f2, '\nSlot-Permeances: Stator=%5.3f,Rotor=%5.3f and Rotor refered to statoi- 
%53f,Lmdss,Lmdsr,Lmddsr); 

fprintf) f2, ’\nSpecific-Slot-Permeance=%5.3f,ssp); 

fprintf) f2, ’\nTotalReactance(X) :Slot+Overhang+Zig-Zag=%5.3f+%5.3f+%5.3f=% 5.3fohms’,Xs,X0,Xz,X); 
fprintf(f2, ’\nShort-Circuit: R=%5.2f,Z=%5.2f and Isc=%5.1f at pl=%5.3f;Isc/IFL=% 5.3f,R,Z,Isc,pfsc,RAT); 

fprintf(f2, '\nTotalLosses(PnL+Pcus+Pcur)=%5.1f+%5.1f+%5.1f=%5.1fWandEfficiency=%5.2f perc',PnL,Pcus,Pcur,Pt,eff); 

fprintf) f2, "vnSlip at FL =%5.3f X Perc'.SFL); 
fprintf(f2, ’VnStarting Tq =%5.3f X FL-Tq'.Tst); 
fprintf(f2, ’\nMax.Output(KW)=%5.1f,Pmax); 
fprintf(f2, ’\nTemp-Rise(deg-C) = %5.1f,Tr); 

fprintf(f2, 

, \nTotalWt(Kg)=Wcus+Wt+Wc+Wri+Wcur+Wcue=%5.1f+%5.1f+%5.1f+%5.1f+%4.1f+%4.1f=%5.1f,Wcus,Wt,Wc,Wri,Wcur,Wcue,Wtot); 
fprintf) f2, ’\nTotal Wt(Kg)=%5.1f and Kg/KW=%5.2f,Wtot,KgPkw); 

fprintf) f2, '\n< End of Output >'); 

end; 


Impact Factor (JCC): 5.6934 


NAAS Rating: 2.45 
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if optim== 1 

if(sn<12) sn=sn+l; 

if eff >=EFFmax EFFmax=eff; end; 

if abs(eff-EFFmax)<=2e-3 M2=sn; end; 

if KgPkw <=minKgPkw minKgPkw=KgPkw;end; 

if abs(KgPkw-minKgPkw)<=0.001 Ml=sn; end; 

if Tr<=minTr minTr=Tr;end; 

if abs(Tr-minTr)<=0.0001 M3=sn;end; 

if IObyl <=minIObyI minIObyI=IObyI;end; 

if abs(10byl-minl0byl)<=0.0001 M4=sn; end; 

if(cond_type== 1 ) 

if(sn==l) fprintf(f2,'\n Rod conductors are used >');fprintf(f2, 

’\n%3d%3d%4.0f%3.0f%4.0f%3d',P,Zs,D,v,L,S);end; 

fprintf(f2, ’\n%2d%4.1fX%3.1f%4.1f,sn,cond_dia,insS,CDSW); 
fprintf(f2, 

'%5.1fX%2.0f%6.3f%5.1f%4.0f%6.2f%5.0f%6.2f%5.2f%5.2f%5.2f%4. If ,Ws,Hs,Btmax,Tr, DO, eff, Wtot, KgPkw, IObyl, PfO,pfsc,SFL); 

fprintf(f2, ’\n \n’); 

end; 

if(cond_type==2) 

if(sn==l) fprintf(f2,’\n Bar conductors are used >');end; 

fprintf(f2, ’\n%2d%3d%3d%4.0f%3.0f%4.0f%3d%4.1fX%3.1f%4.1f,sn,P,Zs,D,v,F,S,Tstrip,Hstrip,CDSW); 

fprintf(f2, 

’%5.1fX%2.0f%6.3f%5.If%4.0f%6.2f%5.0f%6.2f%5.2f%5.2f%5.2f%4.1f,Ws,Hs,Btmax,Tr,DO,eff,Wtot,KgPkw,I0byI,Pf0,pfsc,SFU; 

fprintf(f2, ’\n \n’); 

end; 

end;end;end;end; 

if(sn==12) 

fprintf(f2, ’Design Variant based on Optimization Criteria; ’); 
fprintf(f2, ’\n for Max. Eff. Variant (Sn)= %3d (%5.2f perc) ',M2,EFFmax); 
fprintf(f2,’\n for Min. Kg /KW Variant (Sn)=%3d(%5.2f)', Ml,minKgPkw) ; 
fprintf(f2, ’\n for Min Temp-Rise Variant (Sn)=%3d (%4.2f)’, M3,minTr) ; 
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fprintf(f2, "Vn for Min lOIr ratio Variant (Sn)=%3d (%4.2f)' ,M4,minI0byI) ; 
end; 

end;end; 
fclose (f2) ; 


Impact Factor (JCC): 5.6934 


NAAS Rating: 2.45 



